Kilohertz Quasi-Periodic Oscillation Peak Separation is not Constant 

in the Atoll Source 4U 1608-52 



M. Mendez 1 - 2 , M. van der Klis 1 ' 3 , R. Wijnands 1 , E. C. Ford 1 , J. van Paradijs 1 ' 4 , B. A. Vaughan 5 



ABSTRACT 

We present new Rossi X-ray Timing Explorer observations of the low-mass X-ray 
binary 4U 1608-52 during the decay of its 1998 outburst. We detect by a direct FFT 
method the existence of a second kilohertz quasi-periodic oscillation (kHz QPO) in its 
power density spectrum, previously only seen by means of the sensitivity-enhancing 
'shift and add' technique. This result confirms that 4U 1608-52 is a twin kHz QPO 
source. The frequency separation between these two QPO decreased significantly, from 
325.5 ± 3.4 Hz to 225.3 ± 12.0 Hz, as the frequency of the lower kHz QPO increased 
from 470 Hz to 865 Hz, in contradiction with a simple beat-frequency interpretation. 
This change in the peak separation of the kHz QPOs is closely similar to that previously 
seen in Sco X-l, but takes place at a ten times lower average luminosity. We discuss 
this result within the framework of models that have been proposed for kHz QPO. Beat 
frequency models where the peak separation is identified with the neutron star spin rate, 
as well as the explanations previously proposed to account for the similar behavior of 
the QPOs in Sco X-l, are strongly challenged by this result. 
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1. Introduction 

In the past two years the Rossi X-ray Timing Ex- 
plorer (RXTE) has discovered kilohertz quasi-period- 
ic oscillations (kHz QPOs) in eight een lo w-mass X-ray 
binaries (LMXB; see van der Klis 1998 for a review). 
In almost all cases the power density spectra of these 
sources show twin kHz peaks that move up and down 
in frequency together as a function of mass accretion 
rate, keeping a separation consistent with being con- 
stant (e.g., Strohmayer et al. 1996a ; Wijnands et al. 
19971 ; Ford et al. |l997j ). In some sources a third QPO 
peak has been detected during type-I X-ray bursts, at 
a frequency consistent with the frequency separation 
of the twin peaks (Strohmayer et al. |1996a ) or twice 
that value (Smith, Morgan, & Bradt |1997 ; Wijnands 
& van der Klis |1997| ; Zhang et al. |1996| ; Wijnands 
et al. 1997), indicating a beat-frequency interpreta- 



tion. The third peak, while not strictly constant in 
frequency but varying by up to 0.4%, has been in- 
terpreted as being close to the neutron star spin fre- 
quency (Strohmayer et al. 



1996a) 



Until recently, Sco X-l stood out as the only ex- 
ample where the separation between the two simul- 
taneous kHz QPOs was not constant, but varied by 
40% (van der Klis et al. 1997) as the mass accre- 



tion rate increased, posing a serious challenge to the 
simple beat-frequency interpretation. Using data of 
the 1996 outburst, and a new technique to in crease 
the sensitivity to weak QPOs, Mendez et al. ( 1998| ) 
found the second of the twin peaks in 4U 1608-52, 
and presented evidence for similar variations by 26 % 
of the peak separation in this source, although only 
significant at the 3.5a level. 

In this paper we report on the results from recent 
RXTE observations obtained during the decay of the 
1998 outburst of 4U 1608-52. We observed again two 
QPOs, and show conclusively that the peak separa- 
tion in 4U 1608-52 changes by more than 40%. The 
changes are remarkably similar to those seen in Sco 
X-l. 

2. Observations 

We used the proportional counter array (PCA) on- 
board RXTE to observe 4U 1608-52 between Febru- 
ary 6 and April 25 1998. The observations were trig- 
gered at the peak of the outburst, as measured by the 
RXTE All Sky Monitor (ASM) experiment, and were 
planned to sample the source over the whole decay 
of the outburst. The data set consists of 40 observa- 



tions, each of them with exposures between 600 s and 
19,980 s, and a total observing time of 191 ks. Fig- 
ure ^] shows the ASM light curve during the outburst, 
where we indicate the dates when we carried out the 
pointed PCA observations. 

Beside the two Standard modes that are available 
to all RXTE-PCA observations, we collected data us- 
ing additional modes with high time, and moderate 
energy resolution, covering the nominal 2 — 60 keV 
energy band of the RXTE-PCA. We also included a 
set of burst trigger and catcher modes to record burst 
data at high time and energy resolution. We recorded 
3 type-I X-ray bursts. A timing analysis of these data 
revealed no burst oscillation, with 95 % confidence up- 
per limits to the amplitude of 6 % rms at the peak of 
the bursts, and 20% at the beginning or the end. 

3. Data Analysis and Results 

We used the Standard 2 data and the standard 
PCA background model (Stark et al. |1997| ) for the 
spectral analysis. First we divided the data into 4 
energy bands, 2.0 — 3.5 — 6.4 — 9.7 — 16 keV, and con- 
structed a color-color and hardness-intensity diagram. 
Based on the spectral data, and on the power density 
spectra (described below) in the range 0.01 — 128 Hz, 
we conclude that 4U 1608-52 slowly moved from the 
upper banana, through the lower banana to the island 
state (Hasinger & van der Klis 1989) as the 2 — 60 keV 
count rate decayed from ~ 14, 500 c/s to ~ 260 c/s. 

We also produced background subtracted X-ray 
spectra that were well fitted by means of a model con- 
sisting of a blackbody and a power law with a high- 
energy cut- off. A ssuming a distance of 3.6 kpc (Naka- 
mura et al. 1989 ), during our observations the 2 — 20 
keV luminosity of 4U 1608-52 decayed from 4.8 x 10 37 
erg s _1 at the highest count rate, to 1.2 x 10 36 erg s _1 
at the lowest count rate. Details of the spectral fits, 
color-color diagrams and the analysis of the low fre- 
quency part of the power spectra will be presented 
elsewhere. 

We divided the high-time resolution data into seg- 
ments of 64 s, and calculated a power spectrum 
for each segment preserving the maximum available 
Nyquist frequency. For each of the 40 observations 
we produced an average power spectrum that we 
searched for QPOs at frequencies above >100 Hz. 
None of the observations obtained before March 23 
showed any QPO in this frequency range. The 95 % 
confidence upper limits on the amplitude for a peak 
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with a FWHM of 10, 50, and 100 Hz were 0.8%, 

1.2 %, and 1.6 % rms, respectively on February 6, and 
2.0%, 3.8%, and 4.6% rms, respectively on March 
19. Starting on March 23, when the luminosity of the 
source (2 — 20 keV) had dropped to 1 x 10 37 erg s -1 , 
we detected QPO peaks in 12 observations (Fig. [l]), 
at frequencies between 400 and 1000 Hz. In 5 ob- 
servations, this straight Fourier analysis revealed two 
simultaneous QPOs, while in the remaining 7 obser- 
vations we only detected one QPO peak. In Fig. || we 
show an example of an observation where we detected 
two simultaneous QPOs at ~ 550 Hz (hereafter the 
lower QPO), and at ~ 900 Hz (hereafter the upper 
QPO). 

For those 7 observations where we only detected 
a single peak in the average power spectrum we ap- 
plied the same technique we described in Mendez et 
al. ( 1998 ). We fitted the central frequency of the 
QPO in each segment of 64 s (or longer when occa- 
sionally we had to average some of the 64 s segments 
to detect the QPO), and then shifted the frequency 
scale of each spectrum to a frame of reference where 
the position of the peak was constant in time. Fi- 
nally, we averaged these shifted power spectra. This 
procedure revealed a second QPO in 5 of those 7 ob- 
servations, in 3 occasions at higher frequencies, and 
in 2 observations at lower frequencies than those of 
the peak originally used to align the individual power 
spectra. For the 2 remaining observations (March 31 
and April 5) we did not detect a second peak. The 
95 % confidence upper limits on the amplitude of a 
second (undetected) QPO with a FWHM of 10 Hz, 
50 Hz, and 100 Hz were 3.5%, 4.7%, and 5.3% rms, 
respectively on March 31, and 5.2%, 8.1%, and 10% 
rms, respectively on April 5. 

For the observation of March 25, when the two 
QPOs were more or less equally strong, we measured 
their FWHM in three energy bands, 3.5 — 6.8 — 11.2 — 
14.9 keV; the width of the lower and upper peak were 

8.3 ± 0.4 Hz, and 88 ± 12 Hz, with no significant de- 
pendence on energy. 

In the rest of this paper we only present the results 
of the analysis of those data where, either directly or 
with shifts, we detected two simultaneous kHz QPOs 
in the power spectrum. These are 10 observations 
of the 1998 outburst (Fig. Q, filled circles), plus 2 
observations o f the 1996 outburst (March 3 and 6; 
Mendez et al. 1998| ), yielding a total of 1,336 power 
spectra each of them of 64 s of data. During these 
observations the 2 — 60 keV fractional amplitudes of 



the lower and upper QPO varied from 5.3 % to 9.1 % 
rms, and from 3.3 % to 8.8 % rms, respectively, and 
the FWHM varied from 4.3 Hz to 9.4 Hz, and from 
53 Hz to 173 Hz, respectively. 

We aligned all the 1,336 power spectra using the 
lower QPO as a reference, and we grouped the data 
in 13 sets of ~ 50 to ~ 150 power spectra, such that 
the frequency of the QPO did not vary by more than 
^10 — 20 Hz within each set. Finally, we combined 
these aligned spectra to produce an average power 
spectrum for each set. We fitted all these 13 power 
spectra in the range 256 — 3000 Hz using a func- 
tion consisting of a constant, representing the Poisson 
noise, and two Lorentzians, representing the QPOs. 
The fits were good, with reduced x 2 < 1-1, and the 
significance of both peaks was always > 3a. In Fig- 
ure [| we plot the frequency difference, Av, between 
the upper and the lower QPO as a function of the 
centroid frequency of the lower QPO, uiaw- As ui ow 
increases from ~ 475 Hz to ~ 865 Hz, the separation 
gradually changes from 325.5 ± 3.4 Hz to 225.3 ± 12.0 
Hz, i.e. by 100.2 ± 12.5 Hz. In Figure | we also plot- 
ted Av versus v\ ovl as measured for Sco X— 1 (van der 



Klis et al. 1997). Despite the big difference in lumi- 
nosity, ~ L Edd for Sco X-l and ~ 0.1 X ^Edd for 4U 
1608-52, both sources behave in the same manner. 

In order to get a significant measurement of the 
properties of the upper QPO as a function of the fre- 
quency of the lower QPO, we had to average together 
power spectra from intervals that were very distant 
in time. To check whether this way of averaging the 
data affected our results we selected three contiguous 
intervals for which the frequency of the lower QPO re- 
mained approximately constant, allowing us to mea- 
sure reasonably well both QPOs simultaneously. Dur- 
ing 2624 s, 2624 s, and 3264 s of the observations of 
1998 March 26 and 27, and 1996 March 3 the lower 
QPO remained more or less constant at ~ 600 Hz, 
~ 770 Hz, and ~ 870 Hz, respectively. In the worst 
case (1996 March 3) the significance of the detection 
of the upper peak was 2.7<r. These three intervals 
showed exactly the same trend as seen in Figure |L 
only the error bars were larger. For these three inter- 
vals Av was 298 ± 9 Hz, 277 ± 23 Hz, and 231 ± 13 
Hz, respectively. 

4. Discussion 

We have confirmed the existence of the second kHz 
QPO in 4U 1608-52, and we have conclusively shown 
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that Au, the frequency separation of both peaks, is 
not constant, but changes as a function of the fre- 
quency of the lower kHz QPO. This is the first case 
where this is observed to occur in an atoll source, at 
a luminosity that is far below the Eddington luminos- 
ity. A similar result was obtained previously by van 
der Klis et al. (1997) only for Sco X-l, a Z-source at 
a near-Eddington mass accretion rate. 

Three different models have been proposed to ex- 
plain the kHz QPOs. In the photon bubble oscilla- 
tions model (Klein et al. 1996a, b) the flow of matter 
is funneled onto the polar caps of a highly magnetized 
neutron star accreting at high rates, and the radiation 
pressure becomes locally super-Eddington. The ac- 
cretion column develops radiation-hydrodynamic tur- 
bulence, and energy is transported to the surface 
of the accretion column by photon bubbles, which 
produc e oscillations in the luminosity. Klein et al. 
( 1996b| ) applied this model to explain the 20, 40, and 
60 Hz QPO and the overall shape of the power spec- 
trum observed in the high-luminosity binary pulsar 
GR O J 17 44-28 (Fishman et al. |l995| ; Kouveliotou et 
al. 1996) which contains a neutron star with a high 
magnetic field (B ~ 2.4 x 10 11 G; Cui |1997|) , and to 
the kHz QPO observed in th e luminous Z-source Sco 
X-l (van der Klis et al. 1996 ), which is accre ting m at- 
ter at a high rate (Hasinger & van der Klis 1989). It 
is widely accepted that in the atoll sources the mag- 
netic field is <10 8 ~ 9 G. In 4U 1608-52 during these 
observations the X-ray luminosity was less than 10 % 
of the Eddington luminosity for a 1.4M© neutron star. 
Combining these numbers, and using standard mag- 
netospheric accretion theory, the area of the polecap 
where the matter accreted onto the neutron star was 
between 20 and 60 % of the total stellar surface (see 
Frank, King, & Raine p992L eq. [6.14]). 



It is there- 



fore quite unlikely that the accretion rate was locally 
super-Eddington. This argues, at least for this source, 
against the photon bubble interpretation of the kHz 
QPOs. 



Titarchuk & Mus limoy ( |1997|) , and Titarchuk, La- 
pidus, & Muslimov ( 1998 ) proposed that all the QPOs 
observed in Z and atoll sources might be explained in 
terms of the rotational splitting of oscillation modes of 
a nearly Keplerian accretion disk. One of the predic- 
tions of this model is that the ratio of the frequency of 
the upper to the lower kHz QPO, u upp /u\ ow , is only a 
function of H/R, the ratio of the half-thickness of the 
disk to the radius of the annulus in the disk at which 
the oscillations occur. The ratio u upp /ui ow (Titarchuk 



et al. 1998, eq. [A9] for m = —2 and k = 1) can vary 
from 5 (for H <C R) to 1.37 (for H ~ R). When we 
detected kHz QPOs, the luminosity of 4U 1608-52 
was <10% -^Edd- For these luminosities, and for any 
realistic structure of the disk, one expects H <C R. 



During our observations of 4U 1608-52, 



gradually decreased from 1.69 ± 0.01 at u\ ow = 475 
Hz to a value of 1.26 ± 0.01 at u\ OVJ = 875 Hz, which 
are inconsistent with the predicted value of 5. In Sco 
X-l, despite the large difference in luminosity, and 
therefore large difference of the inferred structure in 
the inner part of the disk, JAi PP /^iow changes in the 
same manner, varying from 1.54 ± 0.01 at 56 5 Hz to 
1.26 ±0.01 at 852 Hz (van der Klis et al |l997|) . More- 



over, in both cases u upp /u\ ovi reaches lower values than 
those allowed by the model. These results show that 
the disk oscillation model as presented cannot account 
for the observed kHz QPOs in these two sources. 

The most widely accepted models for the kHz 
QPOs are beat-frequency models where the upper 
kHz QPO, JAipp, represents the Keplerian frequency 
of the accreting material in orbit around the neu- 
tron star at some preferred radius (van der Klis et al. 
1996), while the lower frequency peak, u\ ovl , is pro- 
duced by the beating of ^ U pp with another frequency, 
us, identified as the spin frequency of the neutron star. 
Strohmayer et al. ( |l996b[ ) proposed that u upp is pro- 
duced at the magnetospheric radius, while Miller et al. 
(1998) proposed that it originates at the sonic radius. 
As 

— ^upp — ^lowj these models predict that the 
frequency difference Au of the twin kHz peaks should 
remain constant, although ^ upp and u\ ow may vary in 
time. This kind of models provides a natural expla- 
nation for the fact that in most of the sources where 
two simultaneous kHz QPOs have been observed Au 
does not change significantly as the frequency of the 
two QPOs varies, while they are also consistent with 
the presence of a third kHz peak that is someti mes de- 
tected during type-I burst s (Str ohmayer et al. 1996a ; 
Smith, Morgan, & Bradt 1997| ) at a frequency near 
u s or 2u s . 

Two explanations have been put forward to ac- 
count for the variable peak separation in Sco X-l, 
both of them related to the inferred near-Eddington 
mass accretion rate in this source. White & Zhang 



( 1997 ) proposed that the photosphere of the neutron 
star might expand by 35 % while conserving its an- 
gular momentum, and therefore slows down as M in- 
creases. Alternatively, at near-Eddington accretion 
the height of the inner disk might increase and the 
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change of t'upp — ^iow might reflect different values of 
^upp a t different heights in the disk (Lamb, private 
communication) . 

However, we know now that in the low luminos- 
ity source 4U 1608-52 Ai/ also changes significantly, 
in contradiction with a beat-frequency interpretation 
involving the neutron star spin. Neither of the above 
explanations can account for this result in 4U 1608- 
52, as it is accreting matter at rates far below the 
Eddington critical rate. One might argue that 4U 
1608-52 and Sco X-l are different, in some respect, 
from the rest of the sources that show kHz QPOs. 
But if this were the case, it would be difficult to ex- 
plain why the properties of the kHz QPOs observed in 
all sources (including both 4U 1608-52 and Sco X-l ) 
are so homogeneous. Moreover, Psaltis et al. (1199^) 
recently showed that for nine other sources where the 
two QPOs have been observed simultaneously, Av is 
consistent both with being constant and with having 
a similar behavior as that observed in Sco X-l and 
4U 1608-52. 

In the beat-frequency model the lower kHz peak is 
expected to be at least as broad as the upper peak, 
as it is generated by a beat between the upper QPO 
and the (coherent) neutron star spin. This is contrary 
to what we observe for 4U 1608-52, where the upper 
peak is intrinsically broader than the lower peak. For 
instance, during 3840 s on 1998 March 25, when the 
frequency of the lower peak varied between 550 Hz 
and 575 Hz, the FWHM of the lower and upper peak 
were 9.3 ± 0.8 Hz, and 95.4 ± 13.8 Hz (2 - 60 keV), 
respectively. A rapidly variable scattering medium 
around the neutron star could cause the broadening, 
and might even be responsible for the frequency shift, 
if it has a (quasi-) periodicity of its own (Lamb, pri- 
vate communication). The QPO at JAi PP , which in the 



sonic model (Miller et al. 199S) is a beaming oscilla- 



tion, would be more sensitive to this than the peak at 
z^iow, which is supposed to be produced by oscillations 
of the luminosity. If this interpretation is correct, it 
also needs to explain why the width of the upper QPO 
does not depend on energy. 

At this point, none of the modified beat-frequency 
models proposed can explain the varying peak sep- 
aration in 4U 1608-52 (and in Sco X-l). In order 
to maintain a beat-frequency model, the change in 
zAipp - fcW by - 40%, both in 4U 1608-52 and Sco 
X-l, requires a similar change in us- The spin of the 
neutron star cannot change by such a large factor on 
such a short timescale. It might still be possible to 



rescue the beat-frequency interpretation, although at 
the expense of its simplicity. Perhaps the idea of a 
layer at the surface of the neutron star that does not 
corotate with the body of the star can still be applied 
to explain these results. The detection of coherent 
oscillations during type-I burst in 4U 1608-52 might 
provide a strong clue to the nature of the kHz QPOs 
in LMXBs, as this would at least make clear at which 
mass accretion level, if any, the kHz peak separation 
does approach the inferred spin rate. 
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Fig. 1. — ASM light curve of the 1998 outburst of 4U 1608-52. The arrows indicate the dates when we performed 
the pointed RXTE observations. Filled circles indicate the observations where we detect, directly or with shifts, 2 
simultaneous kHz QPOs. Open circles and crosses indicate the observations where we detect 1 and no kHz QPOs, 
respectively. 
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Fig. 2. — Power spectrum of a segment of 4224 s starting on UTC 1998 March 25 16:16, for the full energy band of 
the PC A. Both kHz QPOs arc visible. No shift was applied to this power spectrum (see text). On 64 s time scales 
the FWHM of the lower QPO was ~ 5 Hz, but due to the variation of its central frequency during the observation 
it appears broader in the average power spectrum. The upper QPO was broader, and its FWHM on each segment 
of 64 s was ~ 100 Hz. 
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Fig. 3. — The frequency separation between the upper and the lower QPO peak as a function of the freque ncy o f 
the lower QPO for 4U 1608-52 (circles). The crosses are the data obtained for Sco X-l (van der Klis et al. 1997 ). 
The rise of Ai/ at v\ ow ~ 700 Hz in 4U 1608-52 is marginally significant. 
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